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Abstract: We present calculations of spin-relaxation rates 
of alkali-metal atoms due to the spin-axis interaction acting 
in binary collisions between the atoms. We show that for 
the high-temperature conditions of interest here, the spin re- 
laxation rates calculated with classical-path trajectories are 
nearly the same as those calculated with the distorted-wave 
Born approximation. We compare these calculations to recent 
experiments that used magnetic decoupling to isolate spin re- 
laxation due to binary collisions from that due to the for- 
mation of triplet van-der-Waals molecules. The values of the 
spin-axis coupling coefficients deduced from measurements of 
binary collision rates are consistent with those deduced from 
molecular decoupling experiments, and follow a physically 
plausible scaling law for the spin-axis coupling coefficients. 



I. INTRODUCTION 

Spin-exchange optical pumping [^^ of ^He uses spin- 
exchange collisions between '^He atoms and optically 
pumped Rb atoms to produce large quantities of highly 
spin-polarized '^He for a variety of applications, includ- 
ing medical imaging H and spin-polarized targets M. 
The efficiency of polarized '^He production is determined 
by two fundamental rates: the Rb-He spin-exchange rate 
and the Rb spin-relaxation rate. The measured spin- 
exchange rates p] are in fairly good agreement with 
theory 1^. At the high temperatures needed for the 
■^He spin-exchange rates to exceed the wall relaxation 
rates, both Rb-Rb and Rb-He relaxation limit the effi- 
ciency for spin-exchange. Although it had been assumed 
that collisions between alkali atoms rigorously conserve 
the spin polarization, Bhaskar et al. [[7| discovered that 
rapid spin-relaxation in fact occurs in high-density op- 
tically pumped Cs, with a surprisingly large inferred 
spin-relaxation cross section in excess of 1 A^. The cor- 
responding cross section for Rb-Rb relaxation ||5|,p|-[10[ , 
while smaller, still limits the efficiency of '^He production. 
The yet smaller cross section for K-K [g|,|l| relaxation 
suggests that, if technical difficulties are surmounted, K 
may be the optimum partner for spin-exchange with '^He 

A few years ago, we discovered that 1/2 to 2/3 of the 
alkali-alkali relaxation decouples in magnetic fields of a 
few kG pO| , making implausible the interpretation of the 
relaxation exclusively in terms of binary collisions. Re- 
cently, careful magnetic decoupling studies in low pres- 



sure, isotopically pure Rb and Cs samples definitively 
identified the source of the field-dependent relaxation as 
the spin-axis interaction in triplet molecules ||l3| . The re- 
maining alkali-alkali relaxation mechanism at high mag- 
netic field is then presumably from binary collisions. It is 
the purpose of this paper to show that the deduced values 
of the spin-axis interaction from binary collisions are con- 
sistent with the values recently obtained from magnetic 
decoupling studies of triplet molecules |lj] . We hypoth- 
esize a simple scaling law for the second-order spin-orbit 
interaction that explains the relative magnitudes of the 
observed spin-axis interaction strengths. Nevertheless, 
the measured cross sections are in every case at least a 
factor of 10 larger than would be expected from ab ini- 
tio calculations [1^ . Table | contains a summary of the 
existing data on alkali-alkali spin relaxation. 

The spin-axis interaction between two alkali-metal 
atoms is 
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Here the total electron spin of the valence-electron pair 
is S and C is a unit vector lying along the direction of the 
internuclear axis. The coefficient A(i?), a rapidly decreas- 
ing function of interatomic separation R, is currently be- 
lieved to arise from both the direct spin-dipolar coupling 
(averaged over the electron charge distribution) and the 
spin-orbit interaction in second order M]. Accumulat- 



ing evidence, from both high tcrnpcrature [^|14| and 
low temperature experiments il4|-|l(|, suggests that the 
predicted spin-axis coupling (presumably arising almost 
entirely from second-order, spin-orbit interactions) is too 
small in Cs by a factor of 3 to 4, and in Rb by a factor 
of more than 10 [|2|. Put another way, the theoretical 
spin-relaxation cross-sections for Rb are smaller than ex- 
periment by a factor of more than 50. 

This paper carefully documents how the coUisional av- 
eraging of the interaction (m) leads to a spin relaxation 
rate. The somewhat complicated averaging can be done 
exactly within the limitations of the classical path ap- 
proximation, so the origin of the discrepancy cannot be 
due to any inadequacies of the averaging but must lie 
with the potential (|l|), or with the spin-independent in- 
teratomic potential which describes the classical paths, 
or with the neglect of some unknown coUisional relax- 
ation mechanisms other than sudden binary collisions. 
For very cold collisions, where not so many partial waves 



are involved, one could use the distorted-wave Born ap- 
proximation (DWBA) |1^ to account for any limitations 
of the classical-path method, but at the temperatures of 
interest here, as we will show, the results of the classical 
path approximation are within a few percent of those of 
the DWBA. 

Our recent experiments studied the spin-relaxation 
rates as a function of magnetic field for the three alkali 
atoms K, Rb, and Cs, as described in Refs. 0Jl|,|l^,|l|l . 
We assume here that the remaining alkali density depen- 
dent contributions to the relaxation rate at 12 kG arise 
entirely from binary collisions. The data are summarized 
in Table |. 



II. CLASSICAL PATHS 

In a binary collision at temperatures of a few hundred 
Kelvin many partial waves contribute to spin-exchange 
relaxation, so a classical-path treatment should be ade- 
quate. Methods for averaging over all classical-path cross 
sections to obtain spin relaxation rates were given in an 
earlier paper by Walter et al. [gO| (referred to in the 
following as WHW) for collisions between alkali-metal 
atoms and noble-gas atoms. In particular, the anistropic 
magnetic hyperfine interaction between the spin of a no- 
ble gas nucleus and the electron spin of the alkali-metal 
atom, Eq. (3) of WHW has the same tensor symmetry 
as the spin-axis interaction O), and the detailed calcu- 
lations are sufficiently similar to those of WHW that we 
will simply summarize the results here. 

As outlined in WHW and illustrated by Fig. 1 of 
that paper, we can assume that the orbit of the colliding 
pair follows a trajectory governed by the triplet poten- 
tial Vq = ^o(-R)- The time-dependence of the internuclear 
separation R can be found from the equation describing 
conservation of energy and angular momentum, 
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dR / 62 

dt V i?2 Mw^ 



(2) 



Here M is the reduced mass of the pair of colliding 
alkali-metal atoms and w is their relative velocity. For 
the high-temperature experimental conditions of interest 
here, possible changes in direction of the electron spin S 
cause such small changes in the energy or angular mo- 
mentum of the orbital motion that we can neglect them 
and parametrize the orbital energy with the initial rela- 
tive velocity w and the angular momentum by the impact 
parameter b. 

It is convenient to let the time of closest approach of 
the pair be i = 0, so the orbital angle at time t is 



V' = V'W 



cos 



CW-C(O). 



(3) 



The time dependence of the orbital angle can be found 
by numerical quadrature of the equation for conservation 
of angular momentum 



dAp_ 
Itt 



wb 



(4) 



When averaged over a thermal distribution of trajecto- 
ries, collisions in an alkali-metal vapor of atomic number 
density n will cause the mean longitudinal electron spin 
polarization {Sz) of the atoms to relax at the rate 



-{Sz)^-n{vcT)(Sz} 



(5) 



The rate coefficient can be readily calculated by meth- 
ods analogous to those used in WHW for relaxation due 
to the anisotropic magnetic dipole hyperfine interaction. 
The average over all angles of the collisions can be carried 
out analytically and we find in analogy to (33) of WHW, 
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dw p{w) w 



Jo 3 ^^, 



mr, 



(6) 



The probability p(w) dw of finding the magnitude w of 
the relative velocity of the colliding pair between w and 
w + dw is 



p{w) dw = Attw' 



\2TTkBTj 



3/2 



-Mw^/2kBT 



dw, (7) 



where T is the absolute temperature and ks is Boltz- 
mann's constant. The tensor phases accumulated during 
the collision are 



1 

^2m = T 

n 



dt A dom(V') 



(8) 



Here <iom(V') is a Wigner d function, for example, 
dloiip) = (3cos2V' - l)/2 [|l|. In the integrand of 
Eq. (pi) , both the spin-axis coupling coefhcient A and the 
orbital angle ip are functions of time t, obtained from 
numerical integration of Equations (g) and (Q). Since 
ipi—t) = —ip{t), d2mifp) is an even function of t if rn is 
even and odd if m is odd. Also A is an even function of t 
(measured from the time of closest approach) , so (j)2m is 
identically zero iim — ±1. In practice, the rapid decrease 
of A with increasing internuclear separation r means that 
4'2.±2 is so small as to be negligible in practice and so the 
TO = contribution to (p) dominates. 



III. PARTIAL WAVES 

We can also calculate the spin relaxation due to the 
spin-axis interaction (yj) by the distorted-wave Born ap- 
proximation (DWBA), outfined by Newbury et al. [ p^ , 
which we refer to as Newbury in this section. According 
to (47) of Newbury, the rate coefficient corresponding to 
(I) is 



{va) 



dw p{w) w 
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1=0 



{21 + l)Af; (9) 



For the high-temperature conditions of interest here the 
three matrix elements Xu> for spin-flip scattering from 
an initial wave of angular momentum I to a final wave 
of angular momentum I' — 1,1 ±2 are nearly equal. We 
have assumed exact equality to reduce the double sum 
on / and /' to a single sum on Z in (|9|). 

The classical initial relative velocity w is related to the 
asymptotic spatial frequency k of the scattered wave by 



w 



M 

M 



(10) 



In (M) , the integral over impact parameters b which occurs 
in the classical-path expression (|6|) is replaced by a sum 
over partial waves /. The matrix element of the spin- axis 
coupling coefficient between partial waves I and /' is 



An' ^ k gi{r)\{r)gv{r) dr. 
Jo 



(11) 



The wave functions gi (r) are solutions of the Schrodinger 
equation 



dr"^ r 



2 + -TT^o -k']gi^O, (12) 



with the boundary condition as r — > 

gi-^o (13) 

and with the asymptotic boundary condition for kr/l -^ 

oo 



g; -> sin ( fcr - — + 5i 



(14) 



From comparison of (0) with (pi) we conclude that the 
classical-path and partial wave treatment will give prac- 
tically the same answers if 



J^ /2A ^ ^ 



(15) 



In ( |l5| ) we have assumed that / = kb. The initial relative 
energy of the colliding pair is 



21.2 
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(16) 



Numerical solutions to the differential equations (||) and 
(O) are readily obtained. We have confirmed that rela- 
tion (|l5| ) is indeed true, establishing the equivalence of 
the classical-path and partial-wave methods of calculat- 
ing the spin-relaxation rate coefficients. Although the 
two methods give the same results, the classical-path ap- 
proach is much less numerically intensive, requiring only 
the solution of the simple first-order differential equations 
(0) and (0), whereas the partial- wave analysis requires 
solving the second-order equation (E3), with rapidly os- 
cillating solutions. 



IV. COMPARISON OF AB INITIO 
CALCULATIONS TO EXPERIMENT 

In the context of ultracold collisions, Mies et al. p^ ] 
recently published ah initio calculations of A, as the sum 
of two contributions 



A — Ko + Add- 



(17) 



For the heavier alkali-metal atoms, and for small R, they 
found that second-order spin-orbit interactions, repre- 
sented in (y_7|) by Ag^ , were much larger than the term 
Aqq, which describes the direct interaction between the 
magnetic moments of the two valence electrons. Second- 
order contributions analogous to those responsible for 
AgQ are well known from the theoretical literature on 



molecular spectroscopy [E2 
their calculations of A as fo 
results from atomic units): 



|23| . Mies et al. parameterize 
lows (we have converted their 



A = 



4«| 



Ce-P(B-B.s) _ (^^y 



(18) 



Here as is the Bohr radius, /is is the Bohr magneton and 
gs — 2.00232 is the electronic g-factor. The results ( [l^ ) 
of the ab initio calculations are parametrized as follows: 
for Rb, Rs = 5.292 A, C = .001252 and /3 = 1.84 A"!; 
for Cs, Rs = 5.292 A, C = .02249 and (3 = 1.568 A"!. 
The first term in (Hq) represents Ag^ . The second term. 



A. 
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(19) 



represents the magnetic interaction of electrons, taken as 
point particles separated by a distance R. This is an 
excellent approximation at very large R, but should be 
modified, as we show below, at smaller values of R where 
the most important contributions to spin relaxation oc- 
cur. This modification, however, is unlikely to have a 
major impact on the predictions of Ref. |13| ]. 

Table Pshows cross-sections calculated as described in 
Sections^ and [II. For K, we assume only the classical 
spin-dipolar term, because the spin-orbit contribution es- 
timated by scaling from Mies et al. |l^ is negligible. As 
can be seen from Table ffl, the theoretical estimates are 
smaller than experiment by about a factor of 10 for Cs 
and K, and a factor of almost 60 for Rb where the ab 
initio calculations predict that A goes to zero at i? = 5.5 

A. 

In order to describe the effects of spin-relaxation in a 
number of experiments on Cs, where ( |l8| ) predicts relax- 
ation rates that are much too small, the NIST group has 
chosen to multiply Ag^ (the computed contribution to A 
from second-order spin-orbit interactions) by a constant 
value. This assumes that the i?-dependence of the cal- 
culation is correct |14,M. Guided by new experimental 
data, we will discuss similar scaling arguments in Sec- 
tion fvl. 



V. EFFECTS OF WAVEFUNCTION OVERLAP 

ON A„„ 



The expression (^ of A^^ neglects the spatial distri- 
bution of the electron charge. In this section we present 
a simple estimate of the effect of the spatial distribution, 
and find that its neglect cannot be responsible for the 
discrepancy between experiment and theory. 

A simple estimate of the wavefunction of the triplet 
state of an alkali dimcr is 



I*) =^(ri,r2)|x) 

= N[^a{1)vb{2)-Vb{1)va{2)]\x) 



(20) 



where, for example, </3a(1) is a spatial orbital for electron 
1 centered at nucleus A, N a. normalizing factor, and |x} 
is a three-component spinor representing the triplet spin 
state. The matrix element between a final triplet state 
I'^f) and an initial triplet state |^i) of the electronic 
magnetic dipole interaction (at fixed R), is 



"^12 

= ^(x/|Si-(3CC-i)-S2lx. 



(21) 



where 



•^DD — 
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VXri,r2)|2 (22) 



We make the simplifying assumption that ipA and tpB 
can be approximated by the ground-state wavefunction 
of the valence electron of an isolated alkali-metal atom. 
The results are shown in Fig. H^. The principle effect 
of the wavefunction overlap is to reduce the value of A^^ 
as compared to the point dipole approximation, and to 
reduce the predicted ab initio cross section for Rb-Rb to 
5.3 X lO"^*' cm'^, increasing the experiment/theory dis- 
crepency. Use of better electron wave functions than the 
simple form (gy) are unlikely to change the results by the 
orders of magnitude needed to obtain agreement with ex- 
periment. 



VI. SCALING RELATION FOR K 



relaxation rates, we shall assume that Mies et al. have 
correctly identified the second-order spin-orbit interac- 
tion as a major contributor to A, which implies that Ago 
should be proportional to the square of the P\i2^P?,l2 
fine structure splitting Az^, and inversely proportional to 
the valence-electron binding energy i?, as predicted by 
perturbation theory. To obtain the radial dependence of 
AgQ, we shall make the physically plausible assumption 
that it scales as |(/i(r)p, the valence electron density of 
an unperturbed alkali-metal atom at a distance r from 
the nucleus. (The radial dependence of the values of Ago 
calculated by Mies et al. is very nearly that of \(f>{r)\'^.) 
Thus our scaling law is 



A. 



a<^l*MP. 



(23) 



The fine structure splittings At^/c for Cs, Rb, and K are 
554 cni^^, 237.6 cm~^, and 57.7 cm~^, respectively; the 
binding energies E are 3.89 eV, 4.18 eV, and 4.34 eV. 

For the required wavefunctions (j){r), we use the asymp- 
totic expansion of the Coulomb wave function p6[ , 
namely 



ir) = ^ f ^ 



,-r/r: 



(24) 



where the radius r is measured in Bohr radii as- The 
effective principal quantum number n* of the valence 
electron in its ground state is related to the ionization 
energy E (in eV) and the Rydberg i?oo = 13.61 eV by 
n* = {Roc/EY^'^, and the normalization factor N is given 
by TV = [(47r)i/2(„*)3/2p(„*)]-i^ ^l^cre r{x) is the Euler 
gamma function. 



The universal constant il of (23), which has units of 
volume, is deduced from experiments as follows. Re- 
cent experimental studies [Of of the Cs spin relaxation 
in triplet dimers yield a spin-axis coupling |A^''//i| = 
2.79 GHz. This value reflects a thermal average over 
the rovibrational states of the triplet molecules, but 
for simplicity we will take it to be the value of X/h 
at the Cs triplet dimer equilibrium internuclear sepa- 
ration of i? = 12 as, whence, by Eqs. (O) and (|lS|), 



A' 



Cs _ 



2.94 GHz. Thus, using |0|^ = 3.48 x 10" 



R — 12 qb, we obtain 



n = 2880 



(25) 



In this section we show that the spin-axis coupling co- 
efScients deduced from high temperature experiments on 
triplet molecules are consistent with a simple and plausi- 
ble scaling relation that, in turn, accurately predicts the 
relative binary spin-relaxation cross sections. Clearly, the 
contribution A^^ from the direct interaction of the mag- 
netic dipole moments of the electrons is much too small 
to account for observed relaxation rates in the heavier 
alkali-metal atoms. The ab initio calculation of the ad- 
ditional contribution AgQ is also too small. In order to 
test the consistency of the observed molecular and binary 



roughly an atomic volume. 

In Rb, the equilibrium internuclear separation for 
triplet molecules is R = 11.5 as, where |(/)p = 3.23 x 
10-«* ag^, whence Eqs. (H) and (|l|) yield \^^/h = 
294 MHz, in close agreement with the value |A^'^//i| — 
290 MHz deduced from the observation of spin-relaxation 
due to Rb triplet dimers [ [l2| . 

Using the results of Sections O and III and the spin- 
axis coupling coefficients deduced from the above sim- 
ple arguments, we may readily compute the spin relax- 
ation due to the spin-axis interaction in binary collisions 



between alkali-metal atoms. The calculated cross sec- 
tions for K-K, Rb-Rb and Cs-Cs at 400 K, using the 
well-known ab initio interatomic potentials of Krauss 
and Stevens Q to compute the needed trajectories, 
are respectively 4.4 x 10'^° cm^, 1.4 x 10~^^ cm^ and 
1.1 X 10^^^ cm^. As can be seen from Table ||, these val- 
ues agree with experiment to within a factor of 3 for ru- 
bidium and cesium. We note that for potassium it is the 
dipole-dipole interaction that dominates the contribution 
to A, and so the larger discrepancy with experiment for 
K is may be due to the simplified estimate ( p^ ) of A^d ■ 



VII. CONCLUSIONS 

We have shown here that the strength of the spin-axis 
coupling deduced from measurements of binary collisions 
between alkali atoms is consistent with the values de- 
duced from magnetic decoupling of relaxation due to for- 
mation of weakly bound triplet molecules. The deduced 
spin-axis coupling strengths are much larger than pre- 
dicted by ab initio theory. Clearly more reliable theoret- 
ical estimates of the spin-axis coupling are needed. 
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TABLE I. Spin-relaxation cross sections {va)/v (in 10~^* 
cm^) for alkali atoms. Zero-field "cross sections" include con- 
tributions from both molecular formation and binary colli- 
sions. High-field cross sections are assumed to arise solely 
from binary collisions. The theoretical values are calculated, 
at a temperature of 400 K, using estimates of the spin-axis 
coupling strength either from a simple scaling law or from 
Ref. |13]. Also shown is the ratio of experiment to ab initio 
predictions. 
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R(A) 
1. Ab initio calculations of A^, for Rb from Ref. | J13{ , 
and the modified version obtained by spatially averaging the 
magnetic dipole-dipole contribution as described in the text. 
The vertical solid line is at the classical turning point for a 
zero-impact parameter collision at 500 K collision energy. 



